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Very recently, the CDF Collaboration reported the first non-zero measurement of the B s —¥ fi + fj,~ 
branching fraction. The central value of this measurement is more than 5 times of that predicted 
in the Standard Model and, if confirmed, would indicate the existence of new physics. We consider 
the implications of this measurement for the specific case of the Minimal Supersymmetric Standard 
Model (MSSM), and find that it requires large values of tan/3 (J> 30) and favors moderate values for 
the masses of the heavy higgs bosons (mA, mu ~ 300—1200 GeV). We also discuss the implications of 
this measurement for neutralino dark matter, finding that (within the MSSM) regions of parameter 
space in which the lightest neutralino can efficiently annihilate through the pseudoscalar higgs 
resonance (the A-funnel region) are favored. 
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The study of rare decay modes can provide valuable 
probes of physics beyond the Standard Model, not eas- 
ily accessible by other means. Of particular interest are 
the leptonic decays of the B s (sb, sb) and Bd (db, db) 
mesons. In the Standard Model, such decays are dom- 
inated by Z penguin and box diagrams which include 
a top quark loop. As the amplitudes for these processes 
are helicity suppressed and thus proportional to the mass 
of the final state leptons, one might expect B s and Bd 
decays to t + t~ to be most easily measured. Searches in- 
volving tau leptons are very difficult at hadron colliders, 
however, making such decays currently experimentally 
inaccessible. For these and other reasons, the rare de- 
cay modes B s — > /i + /i~ and Bd — > /^ + M~ are among the 
most promising channels with which to constrain or infer 
physics beyond the Standard Model. 

Over the past several years, the CDF [T] and DO [5] col- 
laborations (and more recently, LHCb [3J have reported 
increasingly stringent upper limits on the B s — > /x + /x _ 
branching fraction, steadily moving closer to the value 
predicted in the Standard Model, (3.2 ± 0.2) x 10" 9 g]. 
Very recently, the CDF collaboration reported the first 
measurement of the B s — > fi + fi~ branching fraction in- 
consistent with a value of zero (at the level of 2.8 a). 
Furthermore, CDF's measurement, B(B S —} f.i + fi~) = 
1.8JZo'g x 10~ 8 , favors a central value that is 5-6 times 
larger than predicted in the Standard Model. From this 
result, the CDF collaboration excludes the branching 
fraction predicted by the Standard Model at the 98.1% 
confidence level (C.L.) [5 . 

It has long been appreciated that extensions of the 
Standard Model, including supersymmetry, can lead to 
large enhancements of the B s — > /i + /i~ branching frac- 
tion [51 [7] . In particular, in supersymmetric models with 
large values of tan/3 (the ratio of the vacuum expecta- 
tion values of the two higgs doublets), this branching 
fraction can be as large as ~10-100 times the value pre- 
dicted in the Standard Model [8TU0J. In this letter, we 



explore the implications of CDF's new measurement for 
supersymmetry, focusing for concreteness on the Minimal 
Supersymmetric Standard Model (MSSM). 

The branching fraction for B s — > /i + /i~ can be written 
asQU[T2]: 
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where fg s is the B s decay constant, and to_b s and tb are 
the mass and lifetime of the B s meson, respectively. The 
Wilson coefficients, Cg, C' s , Cp and C' P , describe the 
relevant short distance physics, including any contribu- 
tions from supersymmetric particles. To accommodate 
the large branching fraction measured by CDF, we focus 
on models with large tan/3. In this case, the process of 
B s — > fJ. + n~ is dominated by diagrams with a higgs boson 
(A or H) and a chargino-stop loop, leading to a contribu- 
tion (neglecting QCD corrections) approximately given 
by: 
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where to^ are the masses of the top squarks, mA is the 
mass of the pseudoscalar higgs, Oj is the angle that diag- 
onalizes the stop mass matrix, and Vtb, Vts are elements 
of the Cabibbo-Kobayashi-Maskawa (CKM) matrix. C' P 
and C' s are each suppressed by a factor of m s /rrib and 
thus provide only subdominant contributions. By inspec- 
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FIG. 1: Values of the B s — > fi + fi~ branching fraction over MSSM parameter space as a function of tan/3 (left), uia (center), 
and A^ /j, 2 tan 6 /3m^ 4 Max(p, m st op)~ 4 (right), where m st0 p is the average of the two stop masses. In the lower frames, the 
thermal relic abundance of neutralino dark matter is required to fall within the measured range of the dark matter density (no 
such constraint was imposed in the upper frames). In each frame, we have applied the collider constraints as described in the 
text. The solid horizontal line denotes the central value as measured by CDF, whereas the dotted and dashed lines represent 
the ler and 90% C.L. ranges of this measurement, respectively. 



tion, we see that the dominant supersymmetric contri- In Fig. [I] we show the distribution of values of the 

butions to the B s — > fi + fjT branching fraction scale with B s — ¥ fi + fi~ branching fraction (after imposing the pre- 



fab /3, as well as with m^ 4 from the propagator, /i 2 '• A\ 
from the mass insertions in the sparticle loop, and with 
either /x -4 or my (depending on which mass in the spar- 
ticle loop is heavier). 

To explore this quantitatively, we have used the pro- 
gram Micromegas [13] to calculate the B s — > f-L + n~ 
branching fraction over a large range of the MSSM pa- 
rameter space. In scanning over supersymmetric param- 
eters, we do not assume any particular supersymmetric 
breaking mechanism(s) but instead allow the parameters 
to be chosen independently from one another. We con- 
sider parameters over the following ranges: M\ and slep- 
ton mass parameters up to 2 TeV, /x up to 3 TeV, M 2 and 
rxiA up to 4 TeV, squark and gluino mass parameters up 
to 10 TeV, A t up to 4 TeV (but not in excess of the stop 
masses) and tan/3 up to 70. For all mass parameters, 
we allow both positive and negative values. We require 
that the lightest supersymmetric particle be uncolored 
and electrically neutral, and impose constraints on higgs 
and charged sparticle masses from LEP-II. We also im- 
pose the recent constraints on squark and gluino masses 
from ATLAS (based on 165 pb _1 of data) [14] . and con- 
straints on the tan (3-rriA plane from CMS (based on only 
36 pb _1 ) [15]. Lastly, we impose that the b — > r 4 ^ [TS] 
and b — > sj [17] branching fractions fall within 2er of 
their measured values. At this point, we apply only the 
upper limit on the magnetic moment of the muon, but 
will return to this issue later in this letter. 



viously described constraints) as a function of tan f3 (left), 
vha (center), and A\ p? tan 6 j3 m~ A Max(/i,TO st0 p)~ 4 
(right), where m st0 p is the average of the two stop 
masses. The correlation between this branching fraction 
and tan/3 is particularly striking, essentially requiring 
large to moderate values of this quantity (tan /3 > 30) . 
Although CDF's B s — > [i + [iT branching fraction mea- 
surement also favors moderate or low values of m^, a 
number of acceptable models were found with quite heavy 
rriA- (Note that although models with wia Ss 300 GeV 
and large tan/3 often lead to a large B s —> /i + /i~ branch- 
ing fraction, this combination is inconsistent with the 
constraints from CMS's di-tau searches [T5].) 

In the lower three frames of Fig. [l] we also require 
that the lightest neutralino freezes-out in the early uni- 
verse with a relic abundance within the range inferred by 
WMAP (n CDM h 2 = 0.1120 ±0.0056) |18J. This require- 
ment removes a large fraction of the models found by our 
scan. In particular, the lightest neutralino is predicted 
to be overproduced in the early universe over much of 
supersymmetric parameter space. Regions of parameter 
space which do not suffer from this problem include those 
in which the lightest neutralino efficiently coannihilates 
with a nearly degenerate sparticle, efficiently annihilates 
due to its large couplings (made possible by sizable hig- 
gsino or wino components of its composition), or effi- 
ciently annihilates through the resonant exchange of the 
pseudoscalar higgs boson, A. As the cross section for neu- 
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FIG. 2: Implications of CDF's measurement of B(B S — > 
H + IJ,~) for neutralino dark matter. In the upper frame, we 
compare the mass of the lightest neutralino to the mass of 
the pseudoscalar higgs in models which fall within the 90% 
C.L. range of CDF's measurement, demonstrating that many 
of the models fall near the A-resonance (2m x w itia)- In this 
frame, black (red) points denote models with a bino-like neu- 
tralino (wino or higgsino-like neutralino). In the lower frame, 
we show the spin-independent elastic scattering cross section 
for neutralinos with nucleons, within this same set of models. 
Although the current constraint from XENON-100 [H] ex- 
cludes only a small fraction of these models, next generation 
direct detection experiments will be capable of testing a large 
fraction of the favored parameter space. 



tralino annihilation to bb and t + t~ through A-exchange 
is proportional to tan 2 /3 and \/m\, we expect the A- 
resonance to be highly efficient in many of the models 
favored by CDF's measurement of B[B S — > /j, + /j,~]. In 
in upper frame of Fig. [2j we compare the masses of the 
lightest neutralino and the pseudoscalar higgs in models 
which yield B(B S — > /i + fj,~) within the 90% C.L. region 
as measured by CDF, and which predict a neutralino 
relic abundance consistent with WMAP. We find that 
many of the models in which the lightest neutralino is 
mostly bino-like (black points) fall near the A-resonance 
(2m x sa tua)- Those models with wino-like or higgsino- 
like neutralinos (red points), however, can annihilate effi- 
ciently through chargino and/or neutralino exchange and 
thus need not (and often do not) lie near this contour. In 
the lower frame of Fig. [21 we plot the spin-independent 
elastic scattering cross section of the neutralino with nu- 
cleons, and compare this to the current upper limits from 
the XENON-100 experiment [13]. Although we find that 



direct detection experiments only rule out a small frac- 
tion of those models favored by CDF's measurement of 
B(B S — > [i + ju~), the elastic scattering cross sections pre- 
dicted in many of these models are not far from the cur- 
rent constraints and are expected to fall within the reach 
of next generation experiments. 

Next, we turn out attention to the anomalous mag- 
netic moment of the muon, which has been measured 
to be a° xp = 11659 2080(63) x 10~ n [H|. When com- 
pared to the prediction of the Standard Model, a^ M = 
116591790(65) x 10~ n [21], this measurement consti- 
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To an extent, the supersymmetric contri- 
■> /i + /x~) are correlated. In 
rare decay, large contribu- 



tutes a 3.2cr discrepancy, 

90) x 10" 11 

butions to <5a M and B(B S - 

particular, as with the B s 

tions to Sa^ are found only within the parameter space 

with large tan/3 and with light sparticles. As a concrete 

example, it has been previously shown that within the 

Constrained MSSM (CMSSM) with M 1/2 = 450 GeV, 

M = 350 GeV, A a = 0, and \x > 0, the \a range of <5a M 

translates to a value of B(B S —> fi + /j,~) that is no smaller 

than 1.0 x 10~ 8 [52], illustrating that the measurements 



of Sap and BR(B. 



-> m + m~ 



) each favor overlapping pa- 



rameter space with large values of tan/3. In more general 
terms, for the large values of tan/3 favored by CDF's 
measurement of BR(B S — > /i + /_t _ ), the observed value 
of Sa^ can be accommodated, so long as the masses of 
the contributing slepton and neutralino/chargino are not 
much heavier than a few hundred GeV. 

In summary, we find that to accommodate the large 
B s —} /i + jjT branching fraction reported by the CDF 
collaboration within the context of the Minimal Super- 
symmetric Standard Model (MSSM), one is forced to con- 
sider models with large values of tan /3 and (in most cases) 
modest values of ttia- Relatively large values of A t and 
light stops are also somewhat favored. In much of the 
supersymmetric parameter space favored by this mea- 
surement, neutralino dark matter annihilates efficiently 
through the resonant exchange of the pseudoscalar higgs, 
A, and scatters elastically with nuclei at a rate not far 
below current constraints from direct detection experi- 
ments. 

In the immediate future, the LHC is expected to be- 
come sensitive to the B s — > fi + yT branching fraction 
as reported by CDF. Earlier this year, for example, the 
LHCb collaboration placed an upper limit of 5.6 x 10 -8 
(only a factor of three above CDF's central value) on 
this branching fraction using only 37 pb _1 of data [3]- 
At present, the LHC experiments have each recorded ap- 
proximately 30 times as much data, and should very 
quickly surpass CDF and DO in sensitivity to B s —> 
H + [iT . This will not only provide an opportunity to con- 
firm CDF's measurement, but will also likely enable a 
more precise determination of this quantity. 
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As we were completing this study, Refs. [231 121] ap- 
peared which also discuss CDF's measurement of B s — s- 
fi + fj,~ within the context of supersymmetry. 
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